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1 Préambule 
1.1 Notations 



Pour tout n et m entiers naturels non nuls, et K un corps : 

On appelle M mj „ (K) l'ensemble des matrices à m lignes et à n colones à 
coefficients dans K. On identifie vecteurs colonne (resp ligne) et matrices ne 
possédant qu'une seule colonne (resp ligne). Pour tout M £ M mj „ (K) , on 
désigne par M T sa matrice transposée. 

Pour toute matrice M , de dimension quelconque mais possédant au moins 
deux lignes, (donc y compris les vecteurs colonne), dans un corps quelconque, 
on appellera My la matrice M tronquée de sa première ligne. 

Pour tout K-espace vectoriel E , on note ||.|| : E — > K + la norme euclidienne 
sur E . Sauf dans le cas où E = R , où on la note | . | (valeur absolue) . 

Pour tout espace affine E a associé à un R-espace vectoriel, on note d : E a x 
E a — » M_i_ la distance euclidienne sur E a . 



1.2 Equation diophantienne linéaire 

Soient n un entier naturels non nul. 

Soient a & Mi >n (Z) et b e Z 

On appelle équation diophantienne linéaire l'égalité 

ax = b (1) 
Où x est l'inconnue et est cherchée exclusivement dans M n ,i (Z). 
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On appelle solution de l'équation diophantienne linéaire ax — b tout x G 
M n ,\ (Z) vérifiant (HJ. 

On appelle résolution de l'équation diophantienne linéaire ax — b tout pro- 
cédé permettant de déterminer si l'ensemble de ses solutions est vide ou non. 

1.3 Equation binaire linéaire 

Soient n un entier naturel non nul. 

Soient a G Mi >n (Z) et b G Z 

On appelle équation binaire linéaire l'équation diophantienne linéaire 

ax = b (2) 
Où x est l'inconnue et est cherché exclusivement dans M n ,i ({0, 1})- 

On appelle solution de l'équation binaire linéaire ax — tout x G M n ,i ({0, 1}) 
vérifiant ([2]). 

On appelle résolution de l'équation binaire linéaire ax — b tout procédé 
permettant de déterminer si l'ensemble de ses solutions est vide ou non. 

On appelle système binaire linéaire un système d'équations binaires linéaire. 

On appelle inéquation binaire linéaire l'inégalité 

ax < b 

Où x est l'inconnue et est cherchée exclusivement dans M„ ; i ({0, 1}). 

2 Partition problem 

2.1 Enoncé du problème 

Soit n G N* et soit S = {ai, ...,a n } G 1 n . 

Existe-t-il Si C S et S2 G S tels que : 

- S\ U 5*2 = S et Si (~1 52 = (Si et S2 forment une partition de S ) 

SaeSi a — SaSS 2 a 

Il a été montré NP-complet par Richard M. Karp dans un article de 1972 
(CD- 
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2.2 Réduction du "partition problem" à une équation bi- 
naire linéaire 



Le problème précédent peut encore se formuler comme suit : 
Soit n G N* et soit S = {ai, ...,a n } C Z n . 
Existe-t-il Ni C {1, ...,n} et N 2 C {1, ...,n} tels que : 

- Ni U N 2 = {1, ...,n} et Ni H N 2 = (Ni et N 2 forment une partition de 
{l,...,n}) 

On nomme ce problème (P). 

Montrons que , pour n et 5 tels que définis ci-dessus, (P) est équivalent au 
problème suivant : 

L'équation binaire linéaire (Epp) 

ax = b (3) 

avec a — (2ai, 2a n ) et b — Y^j=i a ii admet-elle une solution ? 

Supposons que (P) admette une solution. Soit Ni et N 2 tels qu'ils forment 
une partition de {1, ...n} et tels que YljeNi a J = Sjgjv 2 a i- 

Xl 

Alors soit x = I ... le M„,i ({0, 1}) tel que Vj G iVijZj = et Vj G 



ATa.Xj = 1 



Nous avons 

n 

ax = 2a>j%j = 2aj = 2 etj 
3=1 je»2 3'eiV 2 

. Et puisque YljeN a i = J2j£N 2 a i ' nous avons que ax = b. C'est-à-dire 
que x est solution de ([3]) (qui admet donc une solution..) 

Réciproquement, si §2$ admet une solution x . Par définition x G Mi <n ({0, 1}) 
et nous pouvons partager en deux les coordonnées Xi,...,x n de x : celles qui 
valent et celles qui valent 1 . 
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Soit Ni C {!,..., n} tel que Vj G Ni.Xj = et N 2 C {l,...,n} tel que 
Vy € N 2 ,Xj = 1 . Ni et N 2 forment clairement une partition de {1, ...,n} et 
nous avons : 

n 

car x est solution de ([3]) . 
Ou encore : 



Soit : 

(P) admet donc une solution. 

Finalement, on a réduit le "partition problem" au problème de la résolution 
d'une équation binaire linéaire. 

De plus, cette réduction est clairement polynomiale. 

3 Résolution d'une équation diophantienne linéaire 
quelconque 

3.1 Enoncé 

Soient n S N*, a — (ai, a n ) G Mi^ n (Z) et b G Z. 

On considère l'équation diophantienne linéaire (E) suivante : 

ax — b 

Pour tout k G {1, n} soit : 
- a le vecteur (ai, ...,a/s) 



j£N 2 i=l 



n 

= 
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- Pk le plus grand diviseur commun de l'ensemble {ai,..., dk}, avec pour 
convention que P\ = a\ . 



On a 



Vfc e {l,...,n}, ^aj-Z = P fe Z 



(4) 



Par conséquent, Vfc S {2, n} , l'équation diophantienne linéaire a k x = Pk 
admet au moins une solution. 



Vfc G {2, n} , soit x 



( \ 



V < k ) 



une solution particulière de l'équation 



diophantienne linéaire a k x = P^. 

Soit enfin la matrice à n lignes et à n — 1 colones 



M 



( a 2 

P> 

P 2 






V o 



a 3 x 



P-2 



r P 3 



Pi Pi 
p 2 Pi 

a 3 x 2 a,4X 2 ° 



Ps 

a. 

p 3 

o 



Pi 

a.4X 3 ° 

Es 
Pi 



P^-2 



P n 
Pn-1 



3.2 Théorème 1 



On conserve les notations énoncées en 3.1 . 
De deux choses l'une : 

- Soit P n ne divise pas b . Dans ce cas (E) n'admet pas de solution. 

- Soit P n divise b . Et dans ce cas, (E) admet une infinité de solutions. De 
plus, dans ce dernier cas, l'ensemble des solutions de (E) est l'ensemble : 

{x p + Mx'\x' e M„_i,i (Z)} 
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Où x 1 



( < \ 



\ < ) 



est une solution particulière de (E) 



Preuve : Supposons que P n ne divise pas b. D'après (jl]) on a Y%=i a j'^ J 
P n Z. Donc (E) n'admet pas de solution. 

On supposera à partir de maintenant que P n divise b. 



Soit alors x p = 



une solution particulière de (E) 



V < J 



Montrons tout d'abord que : Tous les vecteurs de la forme x p + 
Mx', x' <E M n _i )1 (Z) sont solutions de (E). 

Soit x' e M„_i,i (Z) 

On a : a (x p + Mx') = ax p + aMx' = b + aMx' car x p est une solution de 
(E). 

Et aM est le vecteur : 

V P-2 "'" P 2 ' 



^2 



P 2 



P 3 



o„P„_i 



p,l 



Ou encore 

aAf = ( nia2-H2«l 
V P> 



«3 ■ 



-aizf 2 — a2Z 2 =2 +P2 



En utilisant le fait que Vfc G {2, n — 1} , 



Pfc 
x 2 



est solution de l'équa- 



tion diophantienne a k x — Pk , il vient aM = (0, 0). 
D'où aMx 1 = et par conséquent : 

a (x p + Mx') = b 
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Montrons que : Pour chaque solution x de (E) il existe x' G M n _i ; i (Z) 
tel que x — x p + Mx' 



Soit x 
\ x k J 



( X1 \ 

= x 2 
\ x n J 

et x p ' k le vecteur 



une solution de (E). Vfc G , soit x fc le vecteur 

(A \ 

\k ) 



Soit y G Z tel que a" 1 x" 1 = P n -iy ■ Un tel y existe toujours d'après (@|). 
Comme b (x est solution de (E)), nous avons 



P n -iy + a n x n = b 



(5) 



De même, comme x v est solution de (E) , soit y p & Z tel que a" 1 x p ' 71 1 

Pn-lV" 



Nous avons alors : 



De © et © on tire 



P n -iy p + a n Xn = b 



(6) 



Pn-l {y - y P ) = -On (Xn - Xn) 



En divisant par P n : 



(7) 



Comme les entiers relatifs P p 1 et sont premiers entre eux du fait que 
Pn = pgcd (P n -i, On), on peut appliquer le lemme de Gauss à , et en conclure 
que p" 1 divise x n — a;^ . 



Soit alors fc„_i G Z tel que 
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Si y désigne cette fois un entier relatif tel que a n ~ 2 x n ~ 2 = P n -iV (son 
existence est également assurée par (f4|)). Tandis que y p désigne lui un entier 
relatif tel que a n ~ 2 x Pn - 2 = P n -iy v ■ 

En remplaçant à chaque fois dans ax = b, nous avons cette fois les deux 
égalités suivantes : 



P n -2y + a n -\x n -i + a n x n = b 
P n -2y p + fl n -iiLi + a n x p n = b 



Soit : 



P n -2V + On-iXn-i 



Pn-1 



k n -l = b 



P n -2y F + Cin-lX^-i + CL n Xn = b 



D'où l'on tire : 



Pn-2 {y - y P ) + On-l (x n -l ~ + a ™ i^~p^~ kn -^j 



(8) 



Or 



P n -i = a n 1 x Pn - 1 par définition de x^- 1 . En appelant x Pn ~ 1,n 2 le vecteur 
colonne tel que ses coordonnées sont les n — 2 premières de a;*" -1 , on a : 

*n—i — a x -f- a n -\x n _ l 



Soit enfin z € Z tel que a n 2 x Pn - 1 ' n 2 = P n ^ 2 z (z existe d'après (0])). Nous 
avons alors P n -i = P n -2Z + a-n-iXn-'i e t l'égalité ([8]) devient : 



c / P n , / p \ | ( P n -2Z + cin-ix^ 1 , \ n 
"n-2 (2/ - ît) + On-l (a^n-l - + a« I p «n-l I = 

Ou encore : 
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Pn-2 \y-y p + jr z 



( Pi a " X rt-l t 



Enfin, en considérant le fait que P„-i = pgcd(P n -2,a<n-i) nous pouvons 
appliquer de nouveau le lemme de Gauss, cette fois-ci sur l'égalité : 



Pn-2 ( „ , a n 

p — [y - y + -p~ z 



a n-l ( P i a " 3 '"-l u 

- u X n -\ X n _ 1 + K n -l 



Et il vient : p"~^ divise x n -i — x n-i + a " X p~ 1 fc„_i . 
Soit alors fc n -2 £ Z tel que 

_ p a n% n -l u Pn-2 , 

%n—l — x n — l p Kn—1 i p Kn—2 



En répétant le même procédé encore n — 3 fois, on obtient ainsi le vecteur 

/ h \ 

k 2 



colonne x' 



de M n _i t i (Z) tel que : 



( ** \ 

xz 



\ X n ) 



( A \ 



+ 



\ < ) 



( £l 

Pi 






V o 



p 2 

P2 
P 3 





a4a: 3 

P4 



C'est-à-dire 

xt — Xjt + Mtx 
Enfin, en ce qui concerne x\, il est tel que 



a\X\ 



J=2 
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Soit 



a\X\ = b - a,jx? - (q 2 , a n ) x (M T x') 
1=2 

= a\x\ — ((d2) o-n) x Mt) x' 

( L 2 

= aix{ - ( (o 2 , ...,a n ) x ... ] ] x' 

V L„ 

Où L\,...,L n sont les vecteurs lignes de M . 
Nous avons donc : 

/ il \ 

aM = (ai, ...,a„) ... = y^ajLj 
\ L n ) j=i 

Or : nous avons vu dans la première partie de la démonstration que 

aM = (0,...,0) 

n 

- (0,..,0) 



Donc : 



L 2 \ « 
(a 2 ,..,a„) [ ... = y^a 3 £j = -aiLi 



On en déduit que 

/ h 

x\ = x\ — a\L\ 

V k n-l 

Ce qui met fin à la preuve que : 

x = x p + Mx' 

avec x' = | ... | e M n _i,i (Z). 
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3.3 Définition 1 



Toute matrice 



_ 02. 

P'2 
Pl 
Pi 



, P 2 



V o 



a?,x. 2 
Ps 

El 
p 3 




^3 



a^x 



T>3 



Pl 

Pi 



a n x 1 n 1 \ 

pZ 



, où Vfc G 



{2,...,n}, 



est une solution particulière de l'équation diophantienne 



linéaire a k x — Pk est appelée matrice de résolution de (E). 
Remarque : il en existe donc une infinité. 

4 Résolution d'une équation binaire linéaire quel- 
conque 

4.1 Enoncé et notations 

Soient n G N*, a = (ai, ...,a n ) G Mi >n (Z) et b e Z 

On appelle E a l'espace affine canonique sur R de dimension n (c'est-à-dire 
l'espace affine associé à l'espace vectoriel W 1 ). On appelle R le repère affine 
canonique de E a . 

Soit C l'hypercube constitué des points de E a tels leurs coordonnées dans 
R soient toutes dans [0, 1] . On appelle sommet de C tout point de C tel que 
toutes ses coordonnées dans R soient toutes dans {0, 1} . On note S l'ensemble 
des commets de C . 

On considère Y équation binaire linéaire (E) suivante : 

ax = b 



Soit H l'hyperplan affine de E a d'équation dans R : 

aX = b 
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Résoudre (E) revient alors à résoudre le problème : 
H n S est-il vide ou non ? 

S'il est non vide, soit P G H n S . Son vecteur coordonnées X dans R est 
alors solution de (E). 

On désignera donc indifféremment par(E) les deux problèmes. 

De plus , on appelle (Ed) l'équation diophantienne linéaire 

ax = b 



4.2 Définition 2 

Clairement, si (Ed) n'admet pas de solution, alors (E) n'en admet pas non 
plus. 

On appelle cas trivial ce cas de figure. 

4.3 Théorème-définition 2 

( è ^ 

Soit O le point de E a de coodonnées dans R : ... .11 s'agit du centre 

V \ ) 

de l'hypercube C . 

Soit P re f le projeté orthogonal de O sur H . 

On conserve les notations de la section 4 . 

Supposons que l'on ne soit pas dans le cas trivial pour (E) . 

Soient X p une solution particulière et M une matrice de résolution de (Ed) 

Soit alors Rh = (P p , B) , où P p est le point de H de vecteur coordonnées 
X p dans R , et où B est la famille constituée des vecteurs colonnes de M. 

Alors Rh est un repère affine du sous-espace affine H. 

Enfin, le vecteur coordonnées X' re j de P re f dans Rh est égal à : 
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X^ ef = (Mr)- 1 ((X ref ) T - X£) 

Où Mt G M„_i.„_i (Z) est égale à M tronquée de sa première ligne, et les 
deux vecteurs (X re f) T et X^ de M n _i t i (M.) , tous deux obtenus également en 
supprimant la première coordonnée de X re f et de X p respectivement. 



Corrolaire 1 : De plus, puisque l'ensemble des solution de (Eu) est exacte- 
ment l'ensemble : 

{I p + Ml'|l'eM„_ M (Z)} 
(d'après le thérème 1) 
Nous avons : 

Soit P un point de H de vecteur coordonnées X' dans le repère Rh . Il y a 
équivalence entre : 

X' e M„_i,i (Z) et 

X = X v + MX' , le vecteur de ses coordonnées dans R , est solution de 
(Ed) . 



Preuve du théorème-définition 2 : Soit X p une solution particulière de 
(Ed) et M l'une de ses matrices de résolution. 

On appelle P p le point de E a de coordonnées X p dans R . 

Soit B la famille de vecteurs (M\, M„_i) , avec Mj e {1, n — 1} , M,- 
la j ième colonne de M . Alors i? est une base de ~È , direction de H . 

En effet, B compte (n — 1) vecteurs , tous dans ~Ê par construction, et 
son linéairement indépendants (il suffit de regarder la forme de M pour s'en 
convaincre - sous-section 3.3). 

Donc (P P ,B) est un repère affine de H . 

On appelle alors Rjj ce repère. Pour tout point de H , de vecteur coordonnées 
X dans R , on note X' son vecteur coordonnées dans Rh- 
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Soit X re f le vecteur des coordonnées de P re f dans R . 
On a alors 

MX' ref = X ref - XP (9) 

Techniquement, pour calculer X' re j , il nous suffit de résoudre dans R le 
système à n équations et à n — 1 inconnues correspondant à (JÔ]) . 

Dans le cas général, ce type de système n'admet pas toujours de solution. 
Mais comme P re f G H , le système admet pour unique solution X' re f, 

Cela signifie qu'il y a une équation inutile dans ce système. Pour le résoudre, 
il suffit alors de "tronquer" M en lui supprimant sa première ligne. On obtient la 
matrice My de M„_i,„_i (R), triangulaire supérieure, de diagonale strictement 
non nulle (c'est-à-dire dont les coefficients sont. tous non nuls). My est donc 
inversible. Et on appelle respectivement (X re f) T et XZ, les vecteur X re f et X p 
tronqués de leur premier élément. 

Nous avons alors, 

X' ref = (Mt)- 1 ((X ref ) T - X?) 



4.4 Théorème-définition 3 : 

On appelle sommet de C tout point de E a tel que ses coordonnées dans R 
appartiennent toutes à {0, 1} . 

Pour tout repère affine R a d'un espace affine F a , on appelle R^-point de F a 
tout point de F a dont les coordonnées dans R a sont toutes à valeurs dans Z . 

On conserve également les notations, utilisées dans la section 4, pour l'énoncé 
suivant : 

Si l'on n'est pas dans le cas trivial pour (E) , (E) admet une solution si 
et seulement si le iïf^-point le "plus proche" de P re f , au sens de la distance 
euclidienne dans H, appartient à S (c'est-à-dire est un sommet de C). Une 
solution de (JS)est alors donnée par le vecteur des coordonnées dans R de ce 
sommet . 
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Corrolaire-définition 2 : Soit A l'application qui, à un vecteur v à coef- 
ficients dans R de dimension quelconque, lui associe le vecteur constitué des 
entiers les plus proches des coefficients de v , coefficient par coefficient, sauf 
dans le cas où la partie fractionnaire du coefficient vaut \, dans ce cas le choix 
est fait au hasard. 

Le vecteur coordonnées dans Rh du Pf^-point le plus proche de P re f est 
alors A (x' ref ) . 

En effet soit v un vecteur de dimension n par exemple. On note vi....v n ses 
coefficients. Soit Ai, A n les coefficients du vecteur A(v) . 

On a 

n 

(rf(MW)) 2 = E("i- i i) 2 

Supposons qu'il existe un vecteur e à coefficients entiers qui soit plus proche 
de v que A (v) pour la distance euclidienne. 

d(v,e) < d(v,A(v)) 

Alors 

n n 
3 = 1 3 = 1 

Et nécessairement, 3jo G {1, ...n} tel que \vj — e JO | < \vj — Aj \. Ce qui est 
impossible au regard de la définition de A . 

Preuve du théorème-définition 3 : 

Tout d'abord , rappelons le lemme 1 suivant : 

S = {p e E a | P est un P z -point de E a et (d (P, O)) 2 = ^} 

Et : 

VP e E a P lZ -point deP a £ S, (d (P, O)) 2 > J 
Avec d (P, O) la distance euclidienne de P à O dans E a . 
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En effet, soit P s G S et X s = ... le vecteur de ses coordonnées dans R 



Nous avons donc, Vj G {1, ...,n} , Xf G {0, 1}. DoncVj G {1, ...,n} , |Xf - Xf\ = 
\ , où Xf 



77 , OU est la ? ième coordonnée de O dans . D'où 



n n 1 

i=i j=i 



De plus, P s est clairement un P' z -point de E a par définition. 

Réciproquement, soit P G E a de vecteur coordonnées X = .. | dans 
R , tel que : P est un P z -point de E a et d (P, O) 2 = f. 

Supposons que P £ S . Alors, comme P est un P z -point de E a , 3jo G 
{1, .,n} , (X j0 > 2) V (X j0 < -1). Et alors \X jo - Xg\ > 1 + ±. 

Et pour toutes les autres coordonnées de P dans P, j G {1, n} donné, 
nous avons | Xj — Xf | > \ car P est un P z -point de E a (et que par conséquent 
WG{l,...,n},(X J >l)V(X ? <0)). 



Alors, 



3=1 

E te-*?)" + 

je{l,...,n}/{j'o} 

n-1 / l x " 
> — r— + 1 + - 



4 



> 

(ce qui prouve au passage la deuxième partie du lemme 1) 



Ceci est impossible puisque (d(P,0)) 2 = j. 

Donc P G S . Ce qui met fin à la démonstration de la première partie du 
lemme 1. 
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Retour à la preuve : Supposons que l'on ne soit pas dans le cas trivial 
pour la résolution de (E) . 

Alors soit M et X p une matrice de résolution et une solution particulière de 
(E) . P p le point de E a de vecteur coordonnées X p dans R a . Et Rh = (P p , B) 
tel que défini dans le théorème-définition 2, repère affine de H , toujours par le 
théorème-définition 2. 

Soit alors P* le iïf^-point le plus proche de P re j . Comme P* est un R H - 
point , c'est aussi un iî z -point. En effet, soit X'* le vecteur des coordonnées de 
P* dans Rh ■ 

Comme P* est un iîf^-point , X'* G Z n_1 par définition. Alors, d'après le 
théorème 1, X* = X p + MX'* est solution de (E D ) . 

D'où, X* G M nj i (Z) . Et comme X* est le vecteur coordonnées de P* dans 
R d'après la définition de Rh , P* est un 7? z -point de E a . 

Si P* G S , alors (E) admet une solution en tant qu'équation binaire linéaire. 
En effet, on a alors P* G H n S . H n S est donc non vide. De plus , X* est 
clairement une solution de (E) . 

Si P* £ S . Comme P* est un i? z -point de E a , d'après le lemme 1, 
d (O, P*) > j. Comme la droite (O, P re f) est perpendiculaire à H par construc- 
tion de P re f , et que P* E H , les vecteurs P re fP et P re fÔ sont orthogonaux. 

Et il vient : 

(d (O, P*)) 2 = (d (O, P ref )f + (d (P re f,P*)) 2 

Alors 

(d(p ref ,p*)) 2 >^-(d( Pref ,n) 2 

Supposons alors que H D S soit non vide. Soit Q G H n S . On aurait : 

(d(P ref ,0)) 2 = (d(P ref ,0)) 2 + (d(P ref ,Q)) 2 

Car Q G H . 

De plus, comme Q G S , (d(P re f,0)) 2 = j (d'après le lemme 1). 
D'où : 

{d{P ref ,Q)f = - A -{d{P ref ,0)) 2 



17 



Alors 

d(P ref ,P*) >d(P ref ,Q) 

Or , Q est un -Rf^-point de H puisque ses coordonnées dans R vérifient (Eu) 
(corrolaire 1). 

Ceci est impossible, puisqu'alors Q serait un iî^-point strictement plus 
proche de P re f que P* . 

Donc, dans ce cas, HC\S = . C'est-à-dire que (E) n'amdet aucune solution. 

4.5 Algorithme de résolution de (E) : description et com- 
plexité 

Toutes les notations et objets mathématiques introduits en 4.1 sont toujours 
valables. 



4.5.1 Préliminaires 



Calcul de complexité, opérations élémentaires : Nous ne nous intéres- 
serons pas dans les calculs de complexité à des majorations fines du temps de 
calcul ou de la taille en espace mémoire nécessaire, issues de la recherche dans le 
domaine de la complexité, mais aux seules majorations intuitives, qui suffisent 
largement pour ce que l'on veut démontrer. 

Cette façon de faire, dans le but avoué de simplifier les calculs dans ce 
premier article, n'exclut bien évidemment pas de futures approches plus fines 
pour évaluer plus précisément la complexité de l'algorithme. 

De plus, pour ne pas surcharger l'écriture , on ne retiendra à chaque fois que 
le ou les termes dominants dans la majoration. 

Ainsi on considère que, pour deux nombres entiers relatifs z\ et z-i : 

L'addition (ou soustraction) de z\ et zi a une complexité en temps et en 
espace (taille pour stocker le résultat) en 

0(ln (max(|z 1 |,|z 2 |))) 

La comparaison de z\ et z-i a une complexité en temps et en espace en 
0(ln (max(|zi|,|z 2 |))) 
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(dans la mésure où l'on effectue la soustraction et qu'on regarde le signe) 
La multiplication de Z\ par Z2 a une complexité en temps et en espace en 

O (In |*i |+ lu 

Si \z\\ > 1 1 , la division euclidienne de Z\ par z-i , a une complexité en 
temps en 

O (m tel) 

(on cherche le quotient dans {0, \z\\} ou {— \z\\ , 0} par dichotomie, 
puis on calcule le reste par soustraction - cette ) 

Et que le quotient et le reste de la division ainsi calculés ont respectivement 
une complexité en taille de 

O (ln te|) et O (ln te|) 

Ces opérations seront appelées opérations élémentaires. 

De plus, pour les majoration en espace de résultats de calculs, on partira 
souvent du principe suivant : la taille pour stocker un nombre entier relatif 
est logarithmique par rapport à sa valeur. A partir du moment où on est au 
moins en base 2, ce principe est toujours vrai et est à la base de tout calcul de 
complexité. 



Complexité en temps et en espace, complexité de l'algorithme : Si 

un algorithme possède une complexité en temps et en espace différents, alors on 
appellera complexité de l'algorithme le maximum des deux. 



Algorithme d'Euclide : Cet algorithme classique (voir |4]) permet d'obte- 
nir le pgcd de deux nombres entiers relatifs a et fi , \a\ > \/3\ , en calculant 
récursivement les termes de la suite (a n ) ngN , définie par récurence double de la 
manière suivante : 



a = M ,ct x = |/3| 
Vn > 0, a n = qn+2a n +i + a n+2 
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Où, Vn > 0, q-n+2 et a n +2 sont respectivement le quotient et le reste de la 
division euclidienne de a n par a n +i . 

Cette suite est décroissante et positive, et tend vers . Le pgcd de a et /? est 
alors le dernier terme ajv+i non nul de la suite. Pour atteindre ajv+i , il aura 
fallu N itérations (une itération étant le calcul d'un nouveau terme de la suite). 

Le théorème de Lamé stipule que N — O (ln \/3\). Voir |2J. 



Algorithme d'Euclide Etendu (AEE) : Cet autre algorithme classique 
(voir [3]) permet de calculer une solution particulière de l'équation diophan- 
tienne linéaire : 

ax + /3y — pgcd (a, j3) 
Avec a et (3 deux entiers relatifs quelconques, et où x et y sont les inconnues 

e z . 

Et ce en un nombre d'itération TV identique à celui de l'algorithme d'Euclide 
appliqué à a et j3 . 

Soit de nouveau (ct n ) neN la suite par : 

a Q = \a\ ,ai = |/3| 

Vn > 0, a n = q n +2CXn+l + «n+2 

Où, Vn > 0, Çn+2 et a n +2 sont respectivement le quotient et le reste de la 
division euclidienne de a n par a n +i 

Et les suites u = (u„) neN et v = (w„) ngN définies par : 

u = 1, vq = 
ui = 0, t'i = 1 
Vn > 0, u n+2 = u n - q n+2 u n+ i 

V n +2 = Vn - q n +2V n +l 

Alors (ujv+i, uat+i) est une solution particulière de 

ax + f3y = pgcd (a, (3) 



20 



Complexité en taille de l'AEE : Les nombres qui interviennt lors de 
l'AEE sont les termes des suites définies ci-dessus jusqu'à la N ième étape . 

En ce qui concerne les termes de la suite (a„) neN , on a ao — \a\ , et (ce n ) nefil 
est une suite décroissante (et positive), donc les termes suivants sont tous < \a\ 
(et positifs). Leur taille est donc en ln \a\ car ce sont des entiers. 

Pour obtenir un majoration (grossière) des termes u n et v n , on part de la 
majoration suivante : pour tout n > 0, q n < \a\ dans l'AEE (voir le paragraphe 
sur les opérations élémentaires). 

On montre par récurence double sur n > 1 que tous les nombres u n sont 
majorés en valeur absolue par 2 \a\ (et donc par |a| ). 

Comme a est non nul, on a bien \m \ = 1 < 2 \a\° et \u2\ — |ç2| < |a| 1 . 
Puis, supposons que la proposition soit vraie au rang n et n + 1. Alors 
\u n+ 2\ = \u„ - q n +2U n+ i\ < \u n \ + \q n+2 \ \u n+ i\ < 2 1 + 2 \a\ \a\ n < 2 \a\ n+1 
Alors la proposition est vraie aussi au rang n + 2 . Cqfd. 

Nous avons alors que chacun des termes de la suite u à une taille en 
0(Nln\a\) = O (ln \(3\ ln \a\) 

N étant en 0(\a(3) d'après le paragraphe précédent. 
Ou encore, en 

0(ln 2 H) 

car \(3\ < \a\ 

Pour les termes de la suite v la majoration en taille est la même (même 
démonstration). 

Comme il s'agit de suites définies par récurence double, on peut se permettre 
de ne stocker que 2 termes de chacune des 3 suites en même temps, dans le 
déroulement de l'algorithme. 

La taille des termes (çn)„ eN étant en O (ln |a|) , elle est négligeable devant 
la taille des termes de u et v . 

L'espace mémoire nécessaire au stockage des nombres courants dans l'algo- 
rithme AEE est donc de taille 

0(21n 2 H) 
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Complexité en temps de l'AEE : A chaque étape n G {2,...,N + 1} 
(qui consiste à calculer les termes de rang n), on a à faire : 

Une division euclidienne (calcul des termes a n et q n ), de complexité en temps 
en O (In H) 

Deux multiplications et deux additions, de même complexité en temps que 
leurs complexités respectives en espace d'après le paragraphe sur les opérations 
élémentaires, pour le calcul des termes de u et v , soit en O (2 ln 2 \a\) 

Une fois encore, on ne garde que le terme dominant O (21n 2 |a|) 

Comme il y a N étapes = O (ln \/3\). Et comme a > (3 , on peut affirmer que 
le temps global de l'algorithme a une complexité en 



0(21n 3 |a|) 



Conclusion sur la complexité de l'AEE : Elle est en 



0(21n 3 H) 



De plus, le résultat de l'AEE étant les nombre 



UAr+iet v N+ i 



Leur complexité en taille (pour chacun de ces deux nombres) est en 



0(ln 2 |a|) 



4.5.2 Lemme-définition 2 : 



Soit y° = (y°) 



fc€{l,...,7l} 



et x° = 




fc€{l,...,7l} 



les famille de nombre définies 



comme suit : 



y] 



,0 — 



= 1 



Vfc G {1, n — 2} , x n-k+i) es t une solution particulière, déterminée 

par l'AEE, de l'équation diophantienne linéaire de degré 2 : 



Pn — kUn—k ~\~ a n — k+l^n — fe+1 = Pn-k+1 
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Et on pose y® = 1 . Et ( x^x®) , est une solution particulière, fournie par 
l'AEE, de l'équation diophantienne linéaire de degré deux a\X\ + 0,2X2 = P2 . 

Ce qui nous fait au total n—1 appels à l'AEE. 

Alors 

vfce{2,...,n},^^^^xf^^^^^^^•■•'^l ïï «eU°) 

V 1=1 1=2 1=3 l=k-l l=k / 

avec la convention V/' G {1, k} Yl\ = i} yf = 1 , est une solution particulière de 
l'équation diophantienne linéaire : 

k 

ajxj = P k 

Preuve du lemme-définition 2 : Pour k = 2 , on a par définition de (a;?, x 2 ) 
, a\x\ + a 2 X2 = Pi ■ 

Pour k — 3 , a\x\y2 + 0,2x^2 + CL3X3 = P2I/2 + a 3 x 3 — P3 P ar définition de 

(yl4) ■ 

Plus généralement, si on suppose le lemme vrai au rang k < n , au rang k+ 1 
nous avons : 



fe+l k 

3=1 i=j 



Cqfd 



4.5.3 Théorème 4 : 

L'algorithme suivant permet de résoudre l'équation binaire linéaire (E) de 
l'énoncé en 4.1. On conserve les notations introduites dans cette sous-section. 



E a i x j II y ° + a k+i x l+i 



3=1 



E a i x °3 II y° y° + a k+i x k+i 

\3 = 1 1=3 j 

P k y° k + a k+1 x° k+1 

Pk+l 
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Preuve : Elle est contenue dans le paragraphe "description de l'algorithme" 



4.5.4 Description de l'algorithme 

Début de l'algorithme : Si P n ne divise pas b : {Eu) n'admet aucune solu- 
tion (théorème 1), donc {E) non plus . 



Si P n divise b . 

Alors Vfc G {2, n} , on calcule les vecteurs 

( ^ngy? \ 

X 2 IL=2 Vl 



x 1 



V 



4 nf=3 y? 

nk—1 n 
l=fc-iW 



D'après le lemme 2, Vfc G {2, ...,n} , a; fc est solution de a x k = Pk- 

On en déduit que X p = -p-x Pn et M G M„_i „ (Z) constituée du vecteur 

ai 
P 2 

et des vecteurs — X Pk , fc G {2, ...,n— 1} mis en colonne dans 



V o ; 

l'ordre sont respectivement une solution particulière et une matrice de résolution 
de {Eu) (par définition). 

( 1 \ 

Soit Xo = ••■ de M n i (Z). Soit X re f = ta T , où o est le vecteur de 

V \ ) 



l'énoncé et a son transposé, et t = 



C'est-à-dire que X re f est le vecteur coordonnées de P re f , le projeté ortho- 
gonal de O sur H, où O est le point de E a qui a pour coordonnées X dans i? 



On calcule le vecteur 2P ||a|| X re f — 2P a {b — aXo) a T , où P a = Pj i 

plutôt que X re / , pour des raisons de complexité. On calcule aussi X p . 
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D'après le théorème-définition 2, X' re j = (Mt) 1 ((X re f) T — X T ) est le vec- 
teur coordonnées de P re f dans Rh- 

On calcule alors 

2P a \\a\\ 2 X' reî = ({M T )- 1 {{X ref ) T -X p T ))2P a \\a\\ 2 

= (Mt)- 1 (2P a \\a\\ 2 (X ref ) T - 2P a \\a\\ 2 x£) 

= (Mt)- 1 ((2P a \\a\\ 2 X ref ) T -2P a \\a\\ 2 X T ) 

D'après le théorème-définition 3, si le -Rf^-point P* le plus proche de P re f 
est dans S (c'est-à-dire que ses coordonnées dans R sont toutes dans {0,1}), 
alors (E) admet une solution (le vecteur coordonnées de P* dans R) , et dans 
le cas contraire (E) n'admet pas de solution. 

On veut calculer X* ,1e vecteur coordonnées de P* : 

X*'= ^ j = A ( X 'ref) 

Où , A est donnée dans le corrolaire-définition 2 (dans lequel on trouve 
également la justification de ce calcul). 

Pour cela, on effectue la division euclidienne de chaque coordonnée 2P a \\a\\ 2 X' r 
{l,...,n} du vecteur 2P a \\a\\ 2 X' re j: par 2P a \\a\\ 2 , de quotient qj et de reste rj. 

On n'a plus qu'à comparer rj à P a \\a\\ 2 . Si rj < P a \\a\\ 2 , alors l'entier le 
plus proche de -X^e/,j est i s * r j > Pa\\a\\ 2 , alors l'entier le plus proche de 

KefJ est QJ + 1 - 

2 

Dans le cas où rj = P a \\a\\ , les deux entiers qj et qj + 1 sont aussi proches 
de rj l'un que l'autre. 

Pour tout j de {1, ...n} on affecte alors à X* la valeur qj ou + 1 en fonction 
du résultat de la comparaison entre rj et P a \\a\\ 2 (en cas d'égalité le choix est 
fait au hasard et n'a aucune importance). 

Alors on a X* = A (^X' re ^j par construction. 

Et on n'a plus qu'à calculer X* = X' p + MX* pour avoir les coordonnées 
de P* dans R et vérifier que P* e 5 ou non (X*a toutes ses coordonnées dans 
{0, l}ou non. 
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Si oui, on décide : (E) admet une solution , c'est X*. 
Si non, (E) n'admet pas de solution. 



Fin de l'algorithme. 



Notation : Nous appelons cet algorithme : Algorithme Diophantien (AD) 



4.5.5 Complexité de l'AD 

Théorème 5 : La complexité de l'AD est en O (n 4 ln 2 (max a ^)) , où max a ,b = 
max(|6| ,max je{lj ... ;n} \aj\) 

On rappelle que selon les notations, n e N* est le nombre d'inconnues de 
l'équation binaire linéaire à résoudre, b S Z est le membre de droite, et enfin 
di,...a n 6 Z sont les coeffients multiplicatifs affectés aux variables. 



Preuve : On va calculer la taille et le temps de calcul de chaque 
élément calculé dans l'algorithme 

Notations supplémentaires : Soit m a = maxj e { 1; n y \a,j\ 



Les pgcd Pk = pgcd{a\ 1 ...,dfc) , k G {l,...,n} : Les propriétés du pgcd 
faisant qu'on a, ayant posé P\ = ai, Vfc G {1, n} , Pk+i — pgcd (Pk, ak+i) , le 
calcul de tous les Pk s'effectue donc en n — 1 appels à l'AE, de complexité en 
temps : O (lnm a ). 

Taille des Pk : ils sont tous clairement de taille en O (lnm a ) puisque chacun 
d'entre eux divise m a . 
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Vecteurs X Pk , k G {1, ...n] 



Vfc e {2,...n},X I 



X 2 IL=2 Vl 

*fc-i lli=fc-i Wi 



Tous les termes ce!-, j e {1, ...n} et yf, l E {2, ...n} étant issus de l'AEE avec 
comme instance deux nombres majorés en valeur absolue par m a , ils sont tous 
de taille en O (ln 2 m a ), et calculés en un temps en O (m 2 m a ) . L'ensemble de 
ces nombres est le résultat de n — 1 appels à l'AEE, donc le temps total de calcul 
est en O (nln 2 ra a ). L'élément y® = 1 est négligé. 

Pour le calcul des coefficients des X Pk , l'opération la plus complexe est 
clairement celle du terme x\ nlLi Vi ( ou du terme x% YYi = 2 yf)- 

De taille et de temps de calcul en O (inx® + X)^ 1 l n î/; ) > s °it aussi en 
O (nln 2 m a ) d'après le calcul précédent. 

Conclusion : 

Temps total de calcul : en O (2n ln 2 m a ) 

Taille totale en espace à allouer : la taille de chaque X Pk en O (k ln 2 m a ) , 
donc en O (nln 2 m a ) ■ Et comme il y a n X Pk , la taille totale demandée est en 
O (n 2 ln 2 m a ) 



( ~% \ 



Matrice M de résolution de (E) : Elle est constituée du vecteur 



Pi 




V o ) 

des vecteurs — X Pk , k <E {2, ...,n— 1} , auxquels on rajoute le coefficient 
jç^- (entier), mis en colonne dans l'ordre. 

Taille du coefficients maximum en valeur absolue : 
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Le terme le plus grand des X Pk étant en O (n ln 2 m a ) , et comme pour tout 
fc G {2, n — 1}, p *+* est un entier majoré par m a en valeur absolue, nous avons 
que le coefficient le plus grand en valeur absolue sur l'ensemble des vecteurs 
-|r^X Pfc , k e {2,. ..,n-l} est en O (nln 2 m a + lnm a ) . On peut négliger 

lnm devant nln 2 m a pour simplifier l'expression , sans changer la nature de la 
complexité. Et comme les termes et p| sont aussi de taille négligeable par 

rapport à nln 2 m a , on garde O (nln 2 m a ) 
Taille de la matrice : 

Elle est inférieure à sa dimension fois la taille de son plus grand élément en 
valeur absolue. Sa dimension est de n x (n — 1) , d'où une taille en O (n 3 ln 2 m a ) 

Temps de calcul : 

La juxtaposition des colonnes auxquelles on a rajouté un coefficient est une 
opération qui peut être considérée comme nulle 



Vecteur XP = ±-X p " : D' après les calculs précédents , et la complexité 
d'une multiplication, il compte n éléments de taille et de temps de calcul égal 
0(ln|6| + nln 2 rn a ). 

D'où une complexité en temps et en espace de O (nln \b\ + n 2 ln 2 m a ) 

Soit max a ,b = max {\b\ , m a ) , alors la complexité en temps et en espace d'un 
coefficient est en O (n ln 2 max a bj 

La complexité en temps et en espace du vecteur est en O (n 2 ln 2 max a ^) 
Et comme P a est inférieur ou égal à m" -1 on a : 

2P a \\a\\ 2 X p est en O (nln 2 + n 2 lnm a + 2 lnn x m a + n 2 ln 2 max a ,b) donc 
en O (n 2 ln 2 max a ,b) en ne gardant que le terme de plus haut degré (on a m a < 
max a . b ). 
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Vecteur 2P a X ref \\af 



2P a X ref \\a\\ z = P a [ | .. j ||aj| 2 + ( 2b - a ( .. 



Sa complexité en taille et de temps de calcul sont en 
O (n 2 lnm a + n max (2 ln (m a ) , max (ln 2 \b\ , lnm a )) + n lnm a + ln (m a )) 
(voir complexité des opérations élémentaires), et après simplification , en 

O (n 2 ln (max a ,b)) 



Vecteur 2P a \\af X' ref : 

2P a X' ref \\a\\ 2 = (Mt)- 1 (2P a \\a\\ 2 X ref - 2P a \\af X? 



La matrice Mt étant triangulaire (supérieure), on peut résoudre directement 
le système : 

M T (2P a X' ref |H| 2 ) = (2P a X ref \\af - 2P a XP \\af) t 
afin de trouver 2P a \\a\\ 2 X' re ^ . 

X n -i 



Appelons provisoirement X — I ... I le vecteur 2P a \\a\\ 2 X^j et Y 



... ] le vecteur (entier) (2P a X ref \\a\\ 2 - 2P a XP \\af) . 

Yn-X J T 

A chaque étape de calcul, en commençant par la dernière ligne i = n — 1 , 
on calcule 

P ( n ~ 1 
Xi = \Yi- Yl mi à X i 

n \ j=i+i 



avec la convention Y^j= n m i-jXj = . 
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On rappelle que Vi G {1, .., n — 1} , es t un entier. 

Montrons par récurence sur i que chaque Xi ainsi calculé est entier : 

A l'étape i — n — 1 , on a X n _\ = p p " Y n -\ = p p " PaY^-i ■ Comme P n -i 
divise P a = Ylk=i Pk , X n -i est entier, et multiple de Yl^Zi Pk- On n °te X' n _ x 
l'entier relatif tel que X n _\ = ^rife=i Pk) X' n _ 1 

Supposons que pour i G {2, .., n — 2} donné , tous les Xj, j G {i + 1, n — 1} 
, soient des entiers multiples de 111=1 Pk • On note alors Xj l'entier relatif tel 
que Xj = ^rife=i Pk) Xj (hypothèse de récurence) 

Alors tous les Xj ,j G {i + 1, ...,n — ljsont multiples de Ili=i Pk et 

^ = ^ ^ - Ë "'<... - v <j 
= ((n/.) (n>) «-x> u (n*) 

D'où ATi est lui aussi un entier multiple de IIfc=i Pk et alors l'hypothèse de 
récurence est vraie aussi au rang i + 1 

Retour a la complexité : 

D'après les calculs précédent, et d'après la règle sur l'addition et la multi- 
plication, il vient que Y est de complexité en taille et en temps en 

O (n 2 ln 2 max a .b) 

On note également qu'un élément Yi , i G {1, .., n — l}de Y est en 

O (nln 2 maxa.b) 

On rappelle que Mt est de taille (n — l) 2 et que ses coefficients sont tous 
des entiers relatifs majorés en taille par O (nln 2 m a ) ■ 
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On appelle m,ij,(i,j) G {l,...,n— 1} les coefficients de Mr et tom 

m ax(i,j)e{i,...,n-i} 2 

Soit to x = max r | : „ : : \Xj\ 

Pour tout i G {1, n — 1}, soit m l T = maxj e { i; ... ; „_ij ln |X,-| et tot = m\. 
maxjef!,...^!} ln |X,-| 

Pour tout i € {1, ...,n — 1 } , soit 



T,- = lnIX-l = ln 




Nous avons T, < ln 



, d'où 



Tj < ln m a + max (ln \Yi \ , ln (n x tom) + m^ 1 ) 

< liimai,,^ + max (nln 2 max a ,b , ln (n x mu) + rn^ 1 ) 

< Inmaxa^ + max (nln 2 max a ^ , lnn + nln 2 (max a ^) + m 1 ^ 1 ) 



Ce que l'on peut simplifier en négligeant certain termes en 
Ti < nln 2 (maxa.b) + m 1 ^ 1 

D'où : 



m l T = lnmax (ni^t 1 , Tj) 

< max (to^ 1 , nln 2 (max a . b ) + m 1 ^ 1 ) 

< nln 2 (max a ,b) + m 1 ^ 1 
Par récurence immédiate sur i , il vient que 

rriT — m^p < n 2 ln 2 (max a ,b) + m^T 1 
Or, m,T = ln ^p^- + ln \Y n _i\ < lnm a + nln 2 max a ,b 

Et comme X — 2P a X' re ^ ||a|| 2 est un vecteur entier, on en conclut que chacun 

de ses éléments a une taille en O (n 2 ln 2 {max a ,bj) , en négligeant les termes de 
degré inférieur. 

Et 2P a X' re f \\a\\ 2 a une taille en O (n 3 ln 2 (max a ,b)) 
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Vecteur X* : Il est donné par les divisions euclidiennes de chaque coef- 
ficient de 2P a X' re f ||a|| 2 par 2 ||a|| 2 , puis de la comparaison de chaque reste à 




Complexité en temps : Pour chacune des divisions : en O (n 2 ln 2 (max a ,b)) 
d'après le calcul précédent et la complexité en temps d'une division euclidienne. 
Et puisque chaque reste est compris entre et 2 ||a|| 2 la comparaison se fait en 
O (21nm a ), négligeable devant l'autre terme. 

D'où une complexité totale en temps en O (n 3 ln 2 (max a .b)) 

Complexité en taille : pour les coefficients de X* , qui sont les quotients des 
divisions euclidiennes, il sont tous inférieurs en valeur absolue au coeffients de 
2P a X' re j ||a|| 2 . D'où une taille en O (n 2 ln 2 (max a ,b)) pour chacun des termes. 
Et une taille en O (n 3 ln 2 (max a ,b)) pour le vecteur X* en entier. 



Vecteur X* = X p + MX* : M étant une matrice entière de dimen- 
sion n (n — 1) , dont le plus grand de ses coefficients (en valeur absolue) a 
une taille en 0(nln 2 m a ) et sachant que le plus grand des coeffients en va- 
leur absolue de X* est en O (n 2 ln 2 (max a ,b)) , la taille de ses coeffients est 
en O (n (nln 2 m a + n 2 ln 2 (max a ,b))) (il suffit d'écrire le calcul matriciel). On 
retient O (n 3 ln 2 (max a .b)) en négligeant les termes de degré inférieur. Et le 
vecteur en entier est de taille en O (n 4 ln 2 (max a ,b))- H en va de même de la 
complexité en temps puisque additions et multiplications ont même complexité 
en temps et en espace. 

Enfin, X p étant de taille en O (n 2 lnmax a ,b), et sachant qu'on retient uni- 
quement le terme de plus haut degré, la taille du vecteur X* est en O (n 4 ln 2 (max a 
qu'on retient également pour la complexité en temps (pour la même raison). 



Conclusion : Comme n 4 ln 2 (max a .b) est le terme de plus haut degré 
qui majore la complexité en temps de tous les calculs nécessaires pour faire 
fonctionner l'algorithme, on en conclut que l'algorithme est de complexité en 
O (n 4 ln 2 (max a ,b)) 

5 Application au "partition problem" 
5.1 Résolution du "partition problem" 

Théorème 6 : 
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Le "partition problem" , sous sa forme (Epp) étant un cas particulier du pro- 
blème de résolution d'une équation binaire linéaire (E) quelconque, est résoluble 
par l'AD d'après le théorème 5. 



5.2 Complexité 

Théorème 7 : 

Soit n e N* et S = {ai, a„} c Z" 

La complexité du "partition problem" ayant pour instance n et S est en 



( ( ( 7 W\ 

n 4 ln 2 max 2 max (lad) , la,- 1 

l { ' k '))) 



Corrolaire 3 : Puisque n 4 ln 2 ^max ^2 max" =1 (\aj\) , Y^j=x \ a j \j ) est un po- 
lynôme en l'instance du "partition problem" posé précédemment, le "partition 
problem" est dans la classe P . 



Corrolaire 4 : Ceci implique que P = NP . En effet, on a déjà P C NP , 
et tous les problèmes NP sont réductibles polynomialement au "partition pro- 
blem" car il est ./VF-complet (c'est la définition). Or le "partition problem" est 
dans P (théorème 7), donc tous les problèmes NP sont dans P (par réduction 
polynomiale) . 



Preuve du théorème 7 : On prend le "partition problem" sous sa forme 
d'équation binaire linaire (Epp) (voir 12.2]) . Avec les notations précédentes il 
vient : 



- le nombre de variables binaire est n par construction 



le membre de droite dans l'équation : 
absolue par \ a j\ 



En 



est majoré en valeur 



- les coeffients du membre de gauche sont tous majorés en valeur absolue 

par max™ =1 (2 \a,j\) — 2max™ =1 (\aj\) 

Or, il existe un algorithme, l'AD, qui résoud cette équation binaire linéaire 
(théorème 6) , en O ^n 4 ln 2 (max ^2max™ =1 (|aj|) ,Y^j=i l a jl))) ■ cqfd. 
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